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First-principle calculations with dilferent exchange-correlation functionals, including LDA, PBE 
and vdW-DF functional in form of optB88-vdW, have been performed to investigate the electronic 
and elastic properties of two dimensional transition metal dichalcogenides(TMDCs) with the formula 
of MX 2 (M=Mo,W; X=0,S,Se,Te) in both monolayer and bilayer structures. The calculated band 
structures show a direct band gap for monolayer TMDCs at the K point except for M 0 O 2 and WO 2 . 
When the monolayers are stacked into bilayer, the reduced indirect band gaps are found except for 
bilayer WTe 2 , in which direct gap is still present at the K point. The calculated in-plane Young moduli 
are comparable to graphene, which promises the possible application of TMDCs in future flexible and 
stretchable electronic devices. We also evaluated the performance of different functionals including 
LDA, PBE, and optB88-vdW in describing elastic moduli of TMDCs and found that LDA seems to be 
the most qualified method. Moreover, our calculations suggest that the Young moduli for bilayers are 
insensitive to stacking orders and the mechanical coupling between monolayers seems to be negligible. 

Keywords: Transition metal dichalcogenides; Bilayer strnctnres; Elastic properties; 

Electronic strnctnre; First-principles calculation 
PACS: 71.15.Nc; 62.20.de; 62.20.dj; 73.22.-f 

1. Introduction 

The successful isolation of graphene has brought in great revolutions for modern materials 
science The fascinating properties, such as the exceptionally high electron mobility ( ~ 
and unique combination of high modulus (~ 1000 GPa) and tensile strength 
(~ 100 GPa)t^], have made graphene to be one of the most promising candidates for future 
flexible and stretchable electronics. However, the lack of band gap in pristine graphene limits 
its reality application. Due to the presence of a sizeable direct band gap in the visible frequency 
range two-dimensional transition metal dichalcogenides (TMDGs) monolayers have recently 
suggested as an important candidate for electronic and optoelectronic devices Meanwhile, 
TMDGs also exhibit a number of intriguing optical phenomena such as valley-selective circular 
dichroism and the rich physics associated with the valley degree of freedom . Moreover, their 
band gaps are tunable with thickness and different monolayer TMDGs can also be reassembled 
into designer van der Waals heterostructures I®!, which may lead to even richer physics and 
device application. 

Mechanical property of materials such as elastic modulus is of vital importance in any 
practical device applications. It is thus a necessary prerequisite for the integration of TMDGs 
in various devices to obtain a detailed knowledge of mechanical properties of monolayer and 2D 
heterostructures. Using nanoindentation in an atomic force microscope, Bertolazzi et al. have 



investigated in-plane stiffness and breaking strength of suspended monolayer and bilayer M 0 S 2 , 
which indicates that monolayer M 0 S 2 exhibits exceptional mechanical properties comparable to 
stainless steel. In-plane stiffness of monolayer and bilayer M 0 S 2 are found to be 180±60 Nm“^ 
(corresponding to an effective Young’s modulus of 270±100 GPa) and 260 ±70 Nm“^(200± 
60 GPa), respectively. Recently, the elastic modulus of GVD grown monolayer M 0 S 2 and WS 2 
and their Bilayer Heterostructures are also investigated A high 2D elastic moduli of GVD 
monolayer M 0 S 2 and WS 2 (~ 170 Nm“^) are also found. The 2D moduli of their bilayer 
heterostructures are found to be lower than the sum of 2D modulus for each layer, which shows 
a strong mechanical interlayer coupling between the layers 

The importance of mechanical properties in future devices application also intrigues theo¬ 
retical interest. The in-plane stiffness of monolayer M 0 S 2 have been investigated extensively 
[11,12,13,14,15]_ However, the scatted theoretical results from 124 to 146 Nm“^ seems to be lower 
than the experimental value of 170 Nm“^. Recently, the in-plane stiffness of a series of mono- 
layer MX 2 (M=Mo, W; X=S, Se, Te,0) has also been investigated by Kang et al. and Qakr 
et al. Their results suggested that TMDGs with W (O) atom are found to be much stiffer 
in each chalcogenide (metal) group. It should be noticed that the Young modulus of TMDGs 
obtained using PBE functional seems to be lower than experimental measure as shown above, 
it is thus interesting to evaluate the performance of commonly-used functionals to describe the 
mechanical properties of TMDGs. On the other hand, in contrast to the extensively studies for 
monolayers, elastic and mechanical properties of bilayer TMDGs still lack completely. Previous 
study has suggested a transition from direct gap in monolayer to indirect gap in bilayer for 
TMDGs b>i6,i7,i8]^ which implies a strong electronic coupling between different monolayers in 
bilayer TMDGs. Meanwhile, stacking order is also found to affect the electronic band structure 
and absorption spectra It is thus quite interesting to investigate mechanical coupling be¬ 
tween monolayers and to reveal the role of stacking order on elastic property of bilayer TMDs 
theoretically . 

To understand electronic and mechanical coupling of monolayer in bilayer TMDGs, we have 
investigated electronic and elastic properties of eight kinds of TMDGs in both monolayer and 
bilayer structures , in which different functionals are employed for comparison. The detailed 
structure, electronic properties, elastic constants and Young modulus of TMDGs both in mono- 
layer and bilayer obtained using PDA ,PBE and optB88-vdW functional are given and compared 



with available experiment and theory. Moreover, the performance of different functionals in 
describing the elastic properties of TMDCs is also evaluated. 

2. Computational Details 

The present calculations were performed with the Vienna Ab-initio Simulation Package 
(VASP) code within the projector augmented wave (PAW) method. The local density 
approximation (PDA) and Perdew-Burke-Ernzerhof(PBE) are used in our calculation for com¬ 
parison. To take into account van der Waals forces, which are expected to play a crucial 
role in bilayer, the van der Waals density functional in form of optB88-vdW was used. 
Recent study has been evidenced that the optB88-vdW can give much improved results for 
graphene/metal surface Cut off energy and Brillouin zone (BZ) sampling are determined 

after extensive convergence test. The Brillouin zone is simpled by a (18x 18 x 1) Monkhorst- 
Pack grid and a plane-wave basis set with kinetic cutoff energy of 550 eV is applied. The atomic 
positions of those 2D TMDCs have been optimized by using the conjugate-gradient(CG) algo¬ 
rithm. A vaccum of 20 A was set to prevent interactions between periodic images.The elastic 
constants are obtained by strain-energy method. By applying a set of suitable deformations to 
the single unit cell of the 2D TMDCs, we can get the energy-strain curves, and the correspond¬ 
ing elastic constant can be obtained. The detailed relations between applied strain and elastic 
constant are also shown in Ref The strains ranged from -0.05 to 0.05 in steps of 0.01 are 
used in the present calculation. 

3. Results and Discussion 

3.1. Monolayer TMDCs 

3.1.1. Structure properties 

The structure of monolayer TMDCs is shown in Fig. [^(a), in which the M and X atoms 
occupy the hexagonal honeycomb site alternately. Due to the chemical ratio of M:X=1:2, the 
M sublattice layer is sandwiched between two nearby X sublattice layers and forms an X-M-X 
covalently bonded hexagonal quasi-2D lattice as shown in Fig. [^(b), rather than the planar 
construction of the graphene. 



Figure 1: (Color online) Top view (a) and side view (b) of the monolayer TMDCs. (c) and(d) 
are stacking structure for bilayer. Yellow balls and green balls stand for metal and chalcogen 
atoms respectively. 


Table 1: The optimized structural parameters of monolayer TMDCs using LDA,PBE and 
optB88-vdW method, a and d^-x stand for lattice constant and bond length of two X atoms, 
respectively. The PBE results from Reh^^^^ and Reh^^"^] are also listed for comparison. 



LDA 

PBE 

a(A) 

optB88 


Ref. [13] 

LDA 

da;—x(A) 

PBE 

optB88 

M 002 

2.79 

2.83 

2.83 

2.83 

- 

2.44 

2.47 

2.47 

M 0 S 2 

3.12 

3.18 

3.19 

3.18 

3.18 

3.11 

3.13 

3.14 

MoSe2 

3.24 

3.32 

3.32 

3.32 

3.32 

3.31 

3.33 

3.34 

MoTe2 

3.46 

3.55 

3.56 

3.54 

3.55 

3.59 

3.61 

3.62 

WO 2 

2.79 

2.85 

2.83 

2.83 

- 

2.44 

2.48 

2.48 

WS 2 

3.12 

3.18 

3.19 

3.18 

3.18 

3.12 

3.14 

3.15 

WSe2 

3.24 

3.32 

3.32 

3.32 

3.32 

3.33 

3.35 

3.36 

WTe2 

3.47 

3.55 

3.56 

3.55 

3.55 

3.6 

3.63 

3.63 







































The optimized structural parameters of TMDCs with LDA, PBE and optB88-vdW methods 
are listed in Tabletwo previous PBE calculations h3>i4] listed for comparison. We can End 
that the previous calculations are almost identical with our PBE results, which suggests that 
the present results are reliable. It can be clearly seen that the structure seems to be independent 
to metal element but increase monotonously with atomic number of chalcogen. In addition, it 
can be also found that LDA calculation gives a smaller value than PBE, while the results of 
PBE and optB88-vdW are almost the same, which are also more close to bulk crystals. This is 
in agreement with well accepted knowledge that LDA overestimates the binding of materials. 
The similar result of PBE and optB88-vdW indicates van der Waals functional can provide a 
balanced description for systems dominated by not only vdW force but also covalence bonding, 
and van der Waals interaction has few effects on such a monolayer TMDCs. 


3.1.2. Electronic properties 

Now, let’s focus on band structures of monolayer TMDCs. Since the (semi)-local function¬ 
als are known to give similar band structure for materials, we only give the energy band of 
monolayer TMDCs calculated using PBE functional. As shown in Fig. all those eight kinds 
of single layer TMDCs are semiconductors with a distinct band gap, which is quite different 
from the zero-band-gap graphene. For example, the monolayer M 0 S 2 has a direct band gap 
at the K point with a band gap value of 1.62 eV, which is close to the experimental study 
for atomically thin M 0 S 2 layers The direct band gap at the K point can be also found 
in the case of monolayer MoSe 2 and MoTe 2 , and the corresponding band gap is 1.44 eV and 
1.07 eV, respectively. Howerver, for the case of M 0 O 2 , the band gap is indirect with a value 
of 0.98 eV. For monolayer MX 2 (where M stand for Mo and W) with direct energy gap, the 
band gaps decrease when the chalcogenides element becomes heavier i.e. the energy gap order 
is disulhde>diselenide >ditelluride. To analysis the underlying mechanism, we have performed 
Bader charge analysis, which is know to be able to get the charge transfer between different 
atomic species . Taking M 0 X 2 as an example, we can End the total charge of Mo atom in 
monolayer MoX2 is 4.29, 4.52, and 5.06 e, which suggests that the charge transfer from metal 
to chalcogen atom is 1.71, 1.48, 0.94 e. This indicates that the reactivity and binding decrease 
from S to Te, which can also be supported by the fact that the lattice constants decrease with 
the row number of X in Mo dichalcogenides. Clearly, the strong reactivity of chalcogen atom 
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Figure 2: (Color online) The band structures for monolayer TMDCs calculated using PBE. (a)- 
(h) stand for Mo02,MoS2,MoSe2,MoTe2,W02,WS2,WSe2,WTe2 respectively. The Fermi energy 
is marked with light gray dash lines. The band gaps are shown with red arrows. 




















will enhance the hybridization, which will lead to larger band gap of the compounds. In addi¬ 
tion, the band gap for M 0 X 2 (X=Se, S) seems to be smaller than WX 2 , and MoTe 2 and WTe 2 
have quite close gap. However, For monolayer oxide with indirect-band-gap, the band gap of 
WO 2 is much larger than that of M 0 O 2 . The calculated band gaps agree well with other studies 

[13,27] 


3.1.3. Elastic properties 

Elastic property is another important physics quantity focused on the present work. For 2D 
TMDCs, it is possible to measure the elastic constants of these materials from nanoindentation 
experiments using an atomic force microscope However, since the experimental measure¬ 
ment suffers from some uncertainness such as defect, First-Principle calculation can offer an 
important supplement to experiment. 

The elastic constants for monolayer TMDCs calculated with PBE method are reported in 
Table Since the monolayers of 2 H-MX 2 belong to isotropic structures, the linear elastic con¬ 
stant C 22 is equal to Cn. Linear elastic constants of M 0 X 2 and WX 2 decrease monotonically 
with the increase of chalcogen’s atomic number, which means the decrease of the Young mod¬ 
ulus. And for the same chalcogenide atoms, the WX 2 compounds have larger elastic constant 
than M 0 X 2 compounds. The Poisson ratio which is used to describe the lateral deformation is 
calculated by i/ = C' 12 /C'ii . It is remarkable that the Poisson ratios of those 2D TMDCs are 
larger than graphene except for WTe 2 , which means much better lateral contraction property 
than graphene. The calculated Poisson ratio for M 0 S 2 is also close to other studies h2,i4]_ 

For further discussion of elastic property of those monolayer TMDCs, we also evaluate the 
Young modulus E = (Ch — . As shown in Table the Young moduli from different 

methods are close to each other. PBE and optB88-vdW are almost the same, while the results 
of LDA are larger than others. Although there are some numerical differences, the results of 
three methods present the same trends that Young modulus decrease with the atomic number 
of chalcogen and increase from Mo to W compounds. As shown above, charge transfer decreases 
and lattice constant increases in transition metal dichalcogenides MX 2 when X changes from 
S to Te, which suggests that the reactivity and binding decrease from MS 2 to MTe 2 . The 
weakening binding between the metal and the chalcogen atom thus leads to the smaller Young 
modulus. In addition, we can also hnd that the Young modulus for those monolayer TMDCs 


Table 2: The calculated elastic constants of monolayer TMDCs in unit of Nm ^ calculated 
using PBE functional. The data of graphene is also listed for comparison. 



M 0 O 2 

M 0 S 2 

MoSe2 

MoTe2 

WO 2 

WS 2 

WSe2 

WTea 

graphene 

Cii 

232.92 

134.89 

110.7 

84.72 

253.25 

148.47 

122.01 

91.42 

351.6 

C 12 

84.67 

32.03 

25.66 

20.25 

91.32 

31.20 

22.91 

15.40 

64.10 

C 44 

74.12 

51.43 

42.52 

32.23 

80.96 

58.63 

43.03 

38.01 

143.75 

V 

0.36 

0.24 

0.23 

0.24 

0.36 

0.21 

0.19 

0.17 

0.18 


Table 3: The Young modulus of TMDCs in unit of Nm ^ with PBE, LDA and vdW-DF(optB88- 
vdW) method. 



PBE 

LDA 

optB88 

M 0 O 2 

202.1 

221.8 

199.8 

M 0 S 2 

127.3 

141.2 

123.5 

MoSe2 

104.8 

120.9 

105.7 

MoTe2 

79.9 

94.07 

80.7 

WO 2 

220.3 

254.3 

219.7 

WS 2 

141.9 

156.4 

140.7 

WSe2 

114.4 

132.8 

118.2 

WTe2 

88.8 

101.8 

88.4 


are comparable to graphene. For instance, young modulus for M 0 O 2 and WO 2 are almost 
the two-third of the graphene, which is relatively high. M 0 S 2 , MoSe 2 , WS 2 and WSe 2 have 
a value about one-third of the graphene showing excellent elastic property. The calculated 
Young modulus 141 Nm“^ and 156 Nm“^ using LDA for monolayer M 0 S 2 and WS 2 agree well 
with theoretical studies and experiment measurement (180 ± 60 Nm“^ or around 170 
Nm“^ [ 10 ])_ yy-g gjgg LDA calculations gave the closest results to experiment 

among three functionals, which suggest that LDA calculations are more accurate to describe 
the elastic property of monolayer TMDCs. 



3.2. bilayer TMDCs 

3.2.1. Structure properties 

There are different configurations when monolayers are stacked into bilayer structure 
As Fig. [g( c) shown, five stacking orders for the bilayer TMDCs are considered in this paper. 
Since the optB88-vdW can give a very similar in-plane lattice parameter with PBE and the 
interlayer interaction is dominated by vdW forces, we only analysis the results of the bilayer 
TMDCs using optB88-vdW functional. And LDA calculation is also performed for comparison. 
The AAl and AB with almost degenerate energy are found to be the most stable configuration 
for all bilayer TMDCs, which also agrees with earlier calculations. Firstly, we also give the 
structure parameters with five different stacking orders using LDA and optB88-vdW method 
in Table we can find that, the stack of monolayer does not change the monolayer structure 
much. The lattice constants of monolayer and bilayer TMDCs are almost the same for both 
LDA and optB88-vdW method, while the inter-layer distances represented as dm-m varies with 
the stacking order. Both methods predict that the AA and ABl stacking order have a larger 
value of dfn-m than other stacking orders, while the distance for AAl is the smallest, which 
agrees with the stability order of those configurations. The interlayer distances for bilayer 
TMDCs range from 5.2 to 7.8 A according to our calculation. It is clearly shown that the 
structural parameters including in-plane lattice parameter and inter-layer distance obtained 
from LDA are much smaller than that of optB88-vdW for the same bilayer structures. 

3.2.2. Electronic properties 

The band gaps of TMDCs are known to be tunable with thickness jg thus interest 

to investigate the electronic structure of bilayer TMDCs. Since the low-energy AAl and AB 
are found to have similar band structure, we only give the band structures of bilayer TMDCs 
in AAl stacking order in Fig. Unlike the direct band gap of monolayers, indirect band gap 
are found except for bilayer WTe 2 which still has a direct band gap at the K points as shown 
in Fig. 1^ (h). We can easily find that the band gaps have the same variation behaviour with 
the monolayers form MO 2 to MTe 2 . It has to be mentioned that the valence-band maximum 
(VBM) of bilayer WTe 2 at K point is only 0.06 eV larger than the maximum energy for F point, 
which is much smaller than energy difference in other cases. When monolayers are stacked into 
bilayers , the band gaps were found to decrease , which agree well with previous theoretical 
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Figure 3: (Color online) The band structures for monolayer TMDCs calculated using PBE. (a)- 
(h) stand for Mo02,MoS2,MoSe2,MoTe2,W02,WS2,WSe2,WTe2 respectively. The Fermi energy 
is marked with light gray dash lines. The band gaps are shown with red arrows. 















Table 4: The optimized structural parameter of bilayer TMDCs using LDA and optB88-vdW 
method, a stand for lattice constant, d^-m is the distance of the metal atoms between different 


layers. 


Stacking order 


AIB 


AA 


AAl 


AB 


ABl 




a(A) 

dm —m (-^) 

a(A) 

dm —m (-^) 

a(A) 

dm —m (-^) 

a(A) 

dm —m (-^) 

a(A) 

dm —m (-^) 

M 0 O 2 

LDA 

2.79 

4.81 

2.79 

5.30 

2.79 

4.86 

2.79 

4.82 

2.79 

5.28 


optB88 

2.83 

5.28 

2.83 

5.69 

2.83 

5.23 

2.83 

5.20 

2.83 

5.54 

M 0 S 2 

LDA 

3.12 

6.09 

3.12 

6.80 

3.12 

6.01 

3.12 

5.98 

3.12 

6.68 


optB88 

3.18 

6.36 

3.19 

6.75 

3.19 

6.28 

3.19 

6.24 

3.19 

6.71 

MoSe2 

LDA 

3.25 

6.52 

3.25 

7.05 

3.25 

6.33 

3.25 

6.35 

3.25 

7.05 


optB88 

3.33 

6.64 

3.32 

7.14 

3.32 

6.57 

3.33 

6.57 

3.32 

7.09 

MoTe2 

LDA 

3.47 

7.08 

3.46 

7.73 

3.47 

6.85 

3.47 

6.89 

3.47 

7.71 


optB88 

3.56 

7.25 

3.56 

7.78 

3.56 

7.10 

3.56 

7.11 

3.56 

7.75 

WO 2 

LDA 

2.79 

4.88 

2.79 

5.34 

2.79 

4.90 

2.79 

4.87 

2.79 

5.31 


optB88 

2.83 

5.27 

2.83 

5.68 

2.83 

5.22 

2.83 

5.19 

2.83 

5.58 

WS 2 

LDA 

3.13 

6.18 

3.13 

6.81 

3.13 

6.06 

3.13 

6.04 

3.12 

6.69 


optB88 

3.18 

6.37 

3.19 

6.82 

3.19 

6.29 

3.19 

6.26 

3.19 

6.73 

WSe2 

LDA 

3.25 

6.54 

3.25 

7.14 

3.25 

6.39 

3.25 

6.44 

3.25 

7.1 


optB88 

3.33 

6.70 

3.33 

7.15 

3.33 

6.54 

3.33 

6.61 

3.33 

7.14 

WTea 

LDA 

3.47 

7.13 

3.47 

7.76 

3.47 

6.87 

3.47 

6.95 

3.47 

7.72 


optB88 

3.57 

7.27 

3.56 

7.79 

3.56 

7.12 

3.57 

7.12 

3.56 

7.76 


studies 


3.2.3. Elastic properties 

To reveal the evolution of mechanical properties of TMDCs with thickness, we also cal¬ 
culated the elastic constants. Young modulus and Poisson ratio of bilayer TMDCs in Table 

The trend of bilayer Young moduli is very similar to monolayers. A decrease of Young 
modulus is found from MO 2 to MTe 2 and WX 2 are larger than M 0 X 2 . The Young modulus 
of AA and ABl stacking orders are found to be larger than others according to LDA. Mean¬ 
while, we also hnd the differences of Young modulus for different stacking order are within 5%, 
which indicates that the stacking order have little influence on the elastic properties of bilayer. 
Compared with optB88-vdW functional, LDA has a larger Young modulus. Especially, LDA 
method predicts the Young moduli of around 280 Nm“^ for bilayer M 0 S 2 , which is consistent 
with recent experiment 300±13 Table also shows the Poisson ratio of bilayers 

of TMDCs with different conhgurations. The calculated results show that these bilayers have 
close Poisson ratio to the corresponding monolayer. Among different TMDCs bilayer, M 0 S 2 , 








Table 5: Young modulus and Poisson ratio for bilayer TMDCs with different stacking order. 


Young Modulus(Nm 

AIB 

LDA 

optB88 

AA 

LDA 

optB88 

AAl 

LDA 

optB88 

AB 

LDA 

optB88 

ABl 

LDA 

optB88 

M 0 O 2 

439.7 

400.7 

442.3 

399.9 

439.4 

398.9 

438.9 

399.8 

441.2 

405.5 

M 0 S 2 

281.9 

255.7 

285.9 

249.9 

282.0 

252.2 

276.4 

254.0 

283.2 

248.2 

MoSe2 

236.9 

210.4 

241.1 

213.4 

238.2 

212.8 

236.1 

212.2 

240.7 

213.7 

MoTe2 

182.9 

163.0 

189.0 

167.7 

186.6 

162.0 

182.1 

158.6 

186.8 

161.7 

WO 2 

503.4 

461.9 

509.1 

470.8 

504.2 

471.1 

503.9 

468.9 

509.6 

461.9 

WS 2 

314.1 

291.7 

313.2 

279.6 

309.7 

282.2 

309.9 

280.7 

313.7 

282.3 

WSe2 

262.2 

232.4 

265.7 

233.2 

265.1 

227.3 

262.5 

234.3 

266.1 

232.5 

WTe2 

198.8 

172.5 

204.1 

177.9 

203.2 

175.9 

200.1 

172.3 

202.9 

176.9 


Poisson Ratio 

AIB 

LDA 

optB88 

AA 

LDA 

optB88 

AAl 

LDA 

optB88 

AB 

LDA 

optB88 

ABl 

LDA 

optB88 

M 0 O 2 

0.365 

0.364 

0.366 

0.367 

0.365 

0.367 

0.369 

0.37 

0.365 

0.353 

M 0 S 2 

0.212 

0.239 

0.199 

0.253 

0.227 

0.247 

0.218 

0.234 

0.222 

0.256 

MoSe2 

0.201 

0.224 

0.205 

0.232 

0.214 

0.230 

0.211 

0.223 

0.209 

0.229 

MoTe2 

0.202 

0.221 

0.210 

0.218 

0.213 

0.242 

0.217 

0.248 

0.212 

0.245 

WO 2 

0.342 

0.348 

0.337 

0.331 

0.336 

0.329 

0.341 

0.332 

0.340 

0.348 

WS 2 

0.180 

0.194 

0.195 

0.220 

0.201 

0.219 

0.198 

0.213 

0.197 

0.216 

WSe2 

0.173 

0.191 

0.179 

0.200 

0.174 

0.217 

0.174 

0.189 

0.175 

0.204 

WTe2 

0.154 

0.193 

0.155 

0.188 

0.159 

0.199 

0.160 

0.201 

0.158 

0.189 


MoSe 2 and MoTe 2 have similar Poisson ratio with the value of 0.2, which is also independent 
on the stacking order. While the Poisson ratios for bilayer WSe 2 and WTe 2 , are found to be 
much smaller the others. We also notice that optB88-vdW predicts the larger Poisson ratios 
than PDA. 

It should be noticed that the inter-layer interaction is known to be predominated by vdW 
interaction, an accurate description for such a system may need a more advanced treatment 
such as RPA. However, on the basis of the reasonable description for elastic properties of 
TMDCs, PDA method may be seen as an alternative to describe the elastic properties of 
bilayer TMDCs. Recently, Piu et al. [^‘^^have investigated the elastic modulus of CVD grown 
monolayer M 0 S 2 and WS 2 and their Bilayer Heterostructures and found that the 2D moduli 
of their bilayer heterostructures are lower than the sum of 2D modulus of each layer, which 



shows a strong mechanical inter-layer coupling between the layers. However, according to the 
present results, the Young modulus of bilayer structure calculated by both LDA and optB88- 
vdW is almost twice as that of the corresponding monolayer, which shows that the mechanical 
coupling between monolayers seems to be negligible in bilayers. Considered the experimental 
uncertainty, further effort both in experimental and theory is thus suggested to investigate the 
mechanical coupling of monolayers in bilayer. 

4. Conclusion 

We have investigated the band structure and elastic property for eight kinds of two di¬ 
mensional TMDCs in monolayer and different bilayers using different functionals. Our results 
indicate that all those monolayers and bilayers belong to semiconductors and there is a transi¬ 
tion from direct gaps for monolayers (except for M 0 O 2 and WO 2 ), into indirect gaps in bilayers 
(except WTe 2 ). The calculation also shows that the monolayer TMDCs presents excellent elas¬ 
tic property, which are comparable to graphene. Meanwhile, we also hnd that the LDA method 
is an efficient alternative to more advanced method in describing the elastic property of both 
monolayer and bilayer TMDCs. Furthermore, our results also show that the elastic property 
is only influenced slightly by stacking order and the mechanical coupling between monolayers 
seems to be negligible. 
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